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IONOSPHERE TOTAL ELECTRON MEASUREMENTS 
AS EXTRACTED FROM SATELLITE TRACKING DATA 

ABSTRACT 

The Goddard Range and Range Rate System (GRARR) deter- 
mines the range of a spacecraft by measuring the group delay of 
a modulated wave and the range rate by measuring the Doppler 
phase shift of the carrier. Therefore, it has an inherent capa- 
bility to determine the ionospheric dispersive propagation effects 
on the tracking signals which are of equal magnitude but opposite 
sign for the range and range rate measurements. In this paper 
an analytical technique to obtain corrections and the equivalent 
total electron content along the vertical using VHF GRARR data 
and a radial model of the ionosphere is presented. The results 
of computation for a representative pass is given, and an esti- 
mated accuracy is discussed from a statistical view point. Using 
GRARR range, range rate and angle tracking data from Explorer 
41 (also termed Interplanetary Monitoring Platform - G) values of 
total electron content on the order of 10 17 — 10 18 (electrons/ 
meter 2 ) have been obtained. These values are seen to be in good 
agreement with estimates of the total electron content from 
ionospheric profiles using f o F 2 predict data. 
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IONOSPHERE TOTAL ELECTRON MEASUREMENTS 
. AS EXTRACTED FROM SATELLITE TRACKING DATA 

INTRODUCTION 

Demands for accurately locating and predicting the' positions of satellites, in 
determining continental drift and polar motion, earthquake prediction, and geodesy 
in general, impose stringent requirements for correcting spacecraft tracking 
data for the effects of known physical phenomena. 

The interpretation of radio tracking signals will lead to a bias in range and 
range rate calculated from these data if a free-space velocity (i.e., the velocity 
of light in a vacuum) is assumed for the radio signals. This propagation bias 
which is due to the ionospheric plasma is directly proportional to the number of 
electrons along the ray path connecting the satellite and the tracking station. , 

The NASA Goddard Range and Range Rate tracking system is used in two 
frequency bands, namely VHF (~140 MHz) and S-Band (~2 GHz). The VHF sys- 
tem is ideal for both probing the ionosphere as well as verifying ionospheric 
models which can then be applied to the inherently more accurate S-Band system. 
In the Goddard Range arid Range Rate System ranging is accomplished by the use 
of ranging sidetories. The sum of the ranging sidetones phase-modulates the up- 
link signal from the tracking station and the uplink signal in turn phase-modulates 
the downlink signal from the satellite. The phase shift between the ground- 
transmitted sidetones and the satellite transponder-returned sidetones is used 
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as the basic measure of two-way range between the tracking antenna and satellite. 
The range rate measurement consists of the time in seconds required to count a 
fixed number of cycles of the difference frequency between the received and 
transmitted frequency plus a fixed frequency offset. This fixed frequency offset, 
phase locked to the system reference oscillator permits Doppler sign determina- 
tion and is removed during the data processing. The data as shown in Figure 1 
is consequently a measure of the total change in frequency in the round trip of 
the signal from the ground station through the satellite transponder and back 
(reference 1). 

It can be shown from magneto-ionic theory that the ionosphere causes a time 
delay of the sidetone ranging signal and a doppler phase advance of the carrier 
signal. In the GRARR backing system the time delay and phase advance respec- 
tively are given by: 


-2R„ 


T = 

1 


2K 

cf 2 
eq 


(group delay) 


( 1 ) 


2R 

r = 3- + (Doppler phase advance) (2) 

2 c C f 2 

CI eq 


where R T in the absence of calibration bias is the freespace (line of sight) range 
between the tracking station antenna and the satellite, c is the speed of light, 
is the arithmetic mean of the reciprocal uplink and downlink frequency 

' *eq 

squared and K is proportional to the total electron content along the ray path 
connecting the satellite and the ground station. As will be shown, any range 
calibration bias is removed by differentiation. 
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40.3I V 

Also: K = (MKS) 

I y = electron/meter 2 

f - 141 x 10 6 H 7 

c = (2.997925)10® meters/second 

The purpose of this paper is to show the effect of the ionosphere upon the 
range and range rate biases in these two data types and to extract the total 
electron content from GRARR VHF range, range rate, and angle measurements 
using the time derivatives of equations (1) and (2). 

The analysis in this paper is carried out for the Goddard Range and Range 

Rate VHF tracking data taken from Explorer 41 (IMP-G), a spin stabilized sat- 

? . < 

ellite launched 21 June 1969 in eccentric orbit about the earth (see Figure 2) and 
having the following approximate orbital characteristics. 


Semi-major axis 95000 km 

Argument of perigee 203° 

Eccentricity 0.9 < 

Right Ascension of the 
Ascending Node 104° 

Inclination 85° 


The perigee radius of IMP-G is approximately 9500 km (perigee height of 
3122 km), the apogee radius 180,500 km (apogee height of 174,122 km), and the 
period is 3.4 days. 
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ANALYSIS 


Rewriting equations (1) and (2) 

- 2R, 


r = 
1 


2K 

cf 2 

eq 


( 1 ) 


-2R 


T = 
2 


T 2K 
+ 

cf 2 

cq 


(2) 


where R J is the true freespace (line of sight) range between the tracking station 
antenna and the satellite, c is the freespace speed of light, /p is the arithmetic 

eq 

mean of the reciprocal of the uplink and downlink frequencies squared for the 
GRARR VHF transmission, i.e., 


f 2 

eq 


— + — 
f 2 f 2 

l_t dn. 


and K is proportional to the total electron content along the line of sight and is 
given by 

R-, 


K =/ 


T 
/jjd s 


M 



N e e 2 

6 mw 2 

O 


where 


N c = electron density = electrons/ meter 3 
e = electron charge = 1.602 x 10' 19 coulombs 
m = electron mass = 9.11 * 10‘ 31 kilograms 
e 0 = free space dielectric constant = 8.855 x 10' 12 farad/meter 
w = angular frequency = radians/ second 
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Differentiating (1) and (2) and converting' to range rate by (reference (2)) 

- c<t> . - c <fi 

r = - 

<p + 2 Co 2 CO 

where 0 = the angular rate of change of phase (doppler) at the transmit (or side- 


tone) angular frequency. 

For VHF satellites which are spin stabilized, the bias in the range rate k sp 
(reference 3) due to the effect of the spin must be added 



(3) 




where R represents the differentiated range from the sidetone ranging system 
to be distinguished from r which refers to a measurement from the range rate 
system. 

For short data arcs K can be considered a function of the local elevation 

angle E at the approximate height h which is the height of the peak electron 

density of the ionospheric F 2 layer. 

dK = dK.dE (5 , 

dt dE * dt 

Using a spherically symmetric model for the ionosphere in which the gra- 
dient is in the radial direction only we can write (see Figure 3), 


40.3 I v 
sinE 


(6) 
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where I v is the total electron content along the vertical at the point where the 
line of sight intersects the circle of radius a + h 


dK 

dt 


cosE \ . „ (iF 

(40.3) l v 

sin 2 E I dt 


(7) 


dE 


where E and — are related to the ground elevation angle E' and ground eleva- 


HF * 

tion angle rate by 


E = cos -1 


a + h 


cosE 


'] 


dE 


dt L 


1 


I a 


a + h 


cos 2 E 


'j | sinE' 


a + h/\ dt 


dE' 


From (3), (4), and (7) we can write an expression for I as a function of R, 
r , E and E . 


Iv = <f — R _ ksp )( ,/ 2) r |q | .J. 


sin' 


! e\ 


cosE | 


( 8 ) 


This is the basic equation for the determination of I y from range, range rate, 
and angle tracking data used in this analysis. Substituting I y in (3) and (4) we ob- 
tain expressions for the biases in range and range rate due to the ionosphere. 


AR 



40.3I y \ 
sinE | 


(9) 


Af = 



/ cosE 
| sin 2 E 


(40.3) I y - E 


( 10 ) 


For purposes of calculation, the raw range, range rate, and elevation angle 
data was first smoothed by least squares polynomials and then the midpoints of 
the polynomials and differentiated polynomials used for values of R, r , E, and E 


in (8), (9), and (10). 
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It’ should be noted that the ionospheric corrections to the range and range 
rate data are the negatives of (9) and (10). 

Assuming a daytime value for I v of 4 x 10 17 electrons/m 2 and a satellite., 
passing overhead in a 1000 km circular orbit and a VHF frequency of f eq = 141 *10 6 Hz 
we can plot AR and Ar as a function of the ground elevation angle E\ This is 
shown in Figure 4. It can be seen that an appreciable bias results in both data 
types. The VHF GRARR system at a 4 per second sample rate has a range rate 
resolution of better than 0.2 meters/second. With the 20 KHz sidetone the range 
resolution is on the order of 5 to 10 meters depending on the Doppler rate. 


THE PREDICTED UNCERTAINTY OF 1 

Letting AR, Ar,AE, and AE represent errors in R, r ,E, andE, we can write an 
error Al y in I y from (8) . 

AI V = I V (R + AR, r + Ar, E + AE, E + AE)-I V (R, r, E, E) (H) 


or expanding in a Taylor series and retaining only first order partial derivatives 

&I v= (^)^ + (^) Ar+ (^) AE+ (^‘) A ' i (12) 


From (12) we can write an expression for the variance of I v or a 2 as 



where o 2 . , o 2 , a 2 and a 2 - are the variances of R, r ,E, and E respectively, and a 

R f E E 

is the covariance between E and E . It should be noted that each of the partial 
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derivatives is evaluated at some nominal point M , usually taken at the midpoint 
of the data arc. 

EXPERIMENTAL RESULTS 

Goddard Range and Range Rate VHF tracking data taken from Explorer 41 
(IMP-G) on October 6, 1971 at Carnarvon, Australia was analyzed. Using equa- 
tion 8, the total electron content along the vertical was calculated for four data 
arcs. The results can be seen in Table 1 and Figure 5. It can be seen that there 
is excellent agreement between the observed and predicted values. A second 
degree least squares polynomial was used to smooth the range data, and a straight 
line was used to smooth the range rate data. The sampling rate was four per 
second and for the four data arcs approximately 200 points of range and range 
rate data were analyzed. The estimated standard deviation is shown along with 
the predicted value (with an uncertainty of approximately 10% for the local times 
indicated) which was obtained by a NASA/ GSFC ionospheric model developed by 
examining tens of thousands of topside and bottomside profiles and presented at 
the 14th meeting of COSPAR, Seattle, June 1971 by R. B. Bent of DBA Systems, 
Inc., of Melbourne, Florida (reference 4). 
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CONCLUSIONS 

In this paper it has been demonstrated that the Goddard Range and Range 
Rate System has an inherent capability for obtaining the total electron content 
along the ray path connecting the ground tracking station and the satellite. (From 
the determination of the total electron content, corrections can be made to raw 
range and range rate tracking data). These corrections are possible since the 
GRARR system makes independent measurements of group delay (range data) and 
Doppler phase shift (range rate data). This technique is a new one for the mea- 
surement of total electron content in that it combines both group and phase de- 
lays. The correction of raw range and range rate data is important in that more 
precise orbit determinations are possible and estimation of parameters such as 
continental drift, polar motion, etc., by regression analysis techniques can be 
made. 
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Table I 


Experimental Results 


Local 
Time 
h m s 

Number of 
Measurements 
Range and 
Range Rate 

Observed 
I v x 10' 17 
e/m 2 

One Sigma 
x 10 17 
e/m 2 

Predicted 
I v x 10 17 
e/m 2 

7 36 42 

200 

.99 

.13 

1.19 

7 37 32 

200 

1.24 

.19 

1.21 

7 38 22 

200 

.82 . 

.16 

1.21 

7 40 2.5 

204 

.80 

.10 

1.26 


Carnarvon, Australia 
6 October 1971 
Sampling rate 4/second 
VHF Goddard Range 
And Range Rate. 
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Figure 1. Goddard Range and Range Rate Measuring System 



ORBITAL PARAMETERS 

SEMI-MAJOR AXIS 
ARGUMENT OF PERIGEE 
ECCENTRICITY 
RIGHT ASCENSION OF 
ASCENDING NODE 
INCLINATION 
PERIOD 

LAUNCH DATE 


9500 KM 
203° 

0.9 

104° 

85° 

3.4 DAYS 
21 JUNE 1969 


Figure 2. Orbital Characteristics of Explorer 41 (IMP-G) 
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Figure 3. Ionospheric Geometry 
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Figure 4. Spacecraft Tracking Ionosphere Biases 
for 1000 km Altitude Overhead Pass 
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Figure 5. Example of Integrated Content Predicted Versus Observed. 
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